Electron EDM searches based on alkali or alkaline earth bearing molecules 
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We introduce four new molecules — YbRb, YbCs, YbSr + , and YbBa + — that may prove fruitful 
in experimental searches for the electric dipole moment (EDM) of the electron. These molecules 
can, in principle, be prepared at extremely low temperatures by photoassociating ultracold atoms, 
and therefore may present an advantage over molecular beam experiments. Here we discuss prop- 
erties of these molecules and assess the effective electric fields they contribute to an electron EDM 
measurement. 

PACS numbers: 33.15.-e,31.15.A-,11.30.Er 



I. INTRODUCTION 

It is by now widely acknowledged that experiments 
performed on atoms and molecules afford the best chance 
of measuring, or else limiting the possible size of, the 
electric dipole moment (EDM) of the electron p], 0, Q ■ 
A lone valence electron inside an atom with large nuclear 
charge Z will feel an effective electric field enhanced be- 
yond the field that is applied. For example, in atomic 
cesium the apparent field is 132x£ a p P , where £ app is the 
applied electric field [J]. This enhancement arises from 
the field's ability to mix the opposite-parity 6s and 6p 
states of Cs, and a fairly large applied electric field is still 
required. More recent estimates using ab initio meth- 
ods give enhancement factors for Rb and Cs of 25.74 
and 120.53 respectively Experiments on the RbQ 
and CsQ atoms are currently being pursued to take ad- 
vantage of the enhancement factor and long hold times 
achievable in the alkali atom family. 

Thus, to achieve effective electric fields of about 1- 
2 GV/cm one would have to apply fields of ~10- 
40 MV/cm to these atoms. This is not desirable in prac- 
tice. However, one can hold ultracold Rb and Cs atoms 
in traps for a time long compared to the typical molecular 
beam experiment. Therefore, one applies a smaller elec- 
tric field while taking advantage of the long coherence 
times. It is the trade-off between hold time and effec- 
tive electric field size that makes atomic electron EDM 
experiments attractive and worth pursuing. 

Herein lies the rub: one wants a very large applied 
electric field, but with great electric fields come great 
systematic problems. It would therefore be beneficial if 
the s-p mixing could be achieved without the application 
of large laboratory electric fields. This is where molecules 
become interesting, since they can be quite adept at mix- 
ing s- and p-atomic orbitals on their own, without the 
need for large, applied electric fields. The field require- 
ments to polarize the molecule may then be fairly modest. 
This important insight was initially due to Sandars |8j. 

Typically, experiments aimed at measuring the elec- 
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tron EDM in molecules have considered the simplest case 
of molecules with 2 £ symmetry. Then there is only one 
unpaired electronic spin and it has zero orbital angular 
momentum about the internuclear axis. One appealing 
way to achieve fairly robust molecules for this purpose 
is to pair an atom in an s 2 configuration (such as Yb or 
Hg) with fluorine. Fluorine tends to capture an electron 
from any atom it encounters. Thus, when one of the s- 
electrons is captured by F, a 2 £ molecule is formed, as in 
YbFpJ. Moreover, the largely ionic bond is responsible 
for a large effective internal electric field that efficiently 
mixes the atomic s- and p-orbitals of the heavy atom 

In this article we suggest another class of 2 S molecules 
for EDM searches. Our emphasis is to consider molecules 
composed of atoms that can be laser cooled, thus greatly 
reducing the molecular temperature well below that 
available in molecular beam experiments. Specifically, 
we consider the neutral molecules YbRb and YbCs, along 
with the molecular ions YbSr + and YbBa + . In all cases, 
laser cooling of the neutral atoms has been demonstrated 
[HI HH G3 ■ Indeed, recently short-lived ultracold sam- 
ples of metastable YbRb* have been produced via pho- 
toassociation [TI|. This idea is very similar to the pro- 
posed van der Waals molecule CsXe 14 1 . 

These molecules, all of 2 S symmetry in their ground 
states, will have both advantages and disadvantages with 
respect to other molecules that are considered in EDM 
searches. The chief disadvantage is that their elec- 
tronic wave function suggests a sharing, or co-valency, 
of the bond and this somewhat limits the size of the 
effective electric field that they can sustain. We find 
that these fields are typically far smaller (10-100 times) 
than the leading contender — ThO, with an £ c g — 
104 GV/cm[15]. Their main advantage is that they are 
cold enough to be trapped in either optical dipole traps 
or ion traps, leading conceivably to large coherence times 
within which to perform the experiment. In addition, 
these molecules, consisting of two heavy atoms, will have 
smaller rotational constants than the fluorides, hence can 
be polarized in somewhat smaller fields. This is a feature 
that can help reduce systematic effects. Spin-rotational 
effects will also be reduced since they scale inversely with 
the reduced mass of the system; the electron remains 
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fairly well decoupled from the molecular axis in the lower- 
lying rotational levels of the molecules. Another novel 
- though not necessarily beneficial — feature is that 
both atoms contain large- Z nuclei and thus can each con- 
tribute to the effective field. 

In this article we therefore estimate the basic prop- 
erties of these molecules, especially the effective electric 
fields they harbor. Sec. II will discuss the YbRb and 
YbCs molecules and give the details of the potential en- 
ergy surfaces and molecular parameters of the ground 
and excited states of interest. Sec. Ill will present the 
same information for the YbSr + and YbBa + molecular 
ions, as well as the neutrals from which they are derived. 
Sec. IV will detail the effective electric field calculations 
on the 2 S states of interest. Sec. V will present some 
concluding comments. 



II. ALKALI- YTTERBIUM MOLECULES 

Alkali atoms such as Rb and Cs are the workhorses of 
laser cooling, [Iq |. It is the unpaired s-electron that is 
of interest in an electron EDM search since it is what we 
hope to coax into an s-p hybridized atomic state. In the 
introduction we discussed the role of fluorine, which can 
extract an electron from a closed s-shell atom such as Ba 
or Yb. We can accomplish the same feat by pairing an 
alkali atom such as Rb with a closed s-shell atom such 
as Yb because Yb is more electronegative than any alkali 
save H. This is similar to the pairing of Cs with Xe, as 
in [3], except Xe is a closed p-shell atom. Xe serves 
as a polarizing agent on the Cs atom, causing the Cs 
s- and p-orbitals to hybridize. Moreover, the difference 
in electro-negativity is greater between Cs and Yb than 
between Rb and Yb, so that the YbCs molecule might 
be expected to be the more polar of the two. However, 
it is also important to consider the bonding length. A 
larger bond length would work to counteract the charge 
transfer mechanism because the two atoms would not get 
close enough together to allow the atomic wave functions 
to overlap. In order to understand the true details we 
turn to ab initio methods in order to gain understanding 
into the nature of the bonding process. 

For the YbRb and YbCs systems we use the MOLPRO 
suite of codes [TtJ • We start at an inter-atomic distance 
R = 25 ao and move in to R = 4.8 ao. At R—25 oq we 
perform a self consistent field Hartree-Fock (SCF-HF) 
calculation on the closed shell Yb-alkali molecular ion 
system to obtain a starting wave function. We then per- 
form a second SCF-HF on the neutral YbRb molecule 
followed by a multi-configuration self consistent field cal- 
culation (MCSCF) QJ, El- The MCSCF contains the 
ground X 2 T, as well as the excited a 2 n and o 2 £ states 
arising from the alkali np configurations. This is the 
starting configuration for the multi reference configura- 
tion interaction (MRCI) [2(| HH calculation performed 
in each molecular symmetry group. We use an active 
space of {6,3,3,0} with no closed orbitals. We used the 



TABLE I: Molecular parameters for YbRb and YbCs. Bond 
distance r e , dissociation energy (D), vibration constants (a> e 
and WeXe)> rotation constants (B e and D £ ), and vibration- 
rotation mixing constants a e - r e is in atomic units while all 
others are in cm -1 . 



Molecule 


r e 


D 






B e 10" 2 


D e 10" 9 


a e 10~ z 


YbRb X 2 E 


10.22 


194 


12.6 


0.21 


0.99 


24 


1.65 


YbRb a 2 n 


7.44 


3380 


66.4 


0.33 


1.86 


5.9 


1.00 


YbRb & 2 E 


8.52 


942 


36.0 


0.34 


1.42 


8.8 


1.32 


YbCs X 2 E 


10.69 


182 


9.9 


0.13 


0.69 


14 


0.99 


YbCs a 2 n 


7.91 


2710 


50.9 


0.24 


1.27 


3.2 


0.63 


YbCs 6 2 E 


9.41 


490 


20.9 


0.22 


0.90 


6.6 


092 



Stuttgart basis sets and effective core potentials (ECP) 
with core polarization potentials (CPP) for the alkali 
and Yb. We used the ECP68MDFpotential [H for Yb, 
ECP36SDF potential for Rb [HHf and ECP54SDF po- 
tential for Cs[23|]. All of the basis sets associated with 
the ECPs were left uncontracted. The active spaces in 
the calculations were chosen such that the experimental 
excitation energy of the alkali was reproduced to within 
50 cm" 1 . 

The results of the potential calculations are presented 
in Figs. HIl a and b, for the ground state and the two 
lowest-lying dipole- allowed excited states in Hund's cou- 
pling case (a). The inset is an enlargement of the ground 
state showing that, while far shallower than alkali dimcr 
potentials, it is expected to hold bound states. Table U 
gives the relevant molecular parameters such as bond 
length, vibrational constant, dissociation energy, and ro- 
tation constant. The parameters were obtained by fitting 
the ab initio points to a Morse potential. Both YbCs and 
YbRb exhibit a large bonding length greater than 10 a Q 
in the ground X 2 Y, state. Due to the very shallow poten- 
tial and large bond length, these molecules appear to be 
of a van der Waals-type bonding nature, similar to that 
found in CsXe[Ii]. 

In each figure the solid arrow denotes the photoasso- 
ciation step needed to drive the initial atom pair into 
an excited state molecule. For YbRb, this step has 
already been demonstrated experimentally [La ]. Also 
shown schematically as a dashed arrow is the transition 
necessary to create molecules in the ground electronic, 
vibrational and rotational state. Depending on which 
levels can ultimately be accessed in the first photoassoci- 
ation step, this second transition may not be possible in a 
single step. The quest to produce ultracold, ground state 
polar molecules from ultracold atoms is currently a ma- 
jor research area, which has recently seen the production 
of molecular ground states [25l. [2r3|. 

The ab initio calculations show that while the dom- 
inant contribution to the unpaired electron is from the 
Rb 5s (Cs 6s) orbital, there are appreciable contribu- 
tions from the Rb 5p (Cs 6p) as well as the Yb 6s and 6p 
atomic configurations. A secondary contribution comes 
from a similar arrangement except that the Yb 6p atomic 
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FIG. 1: (Color online) Potential energy surfaces for the YbRb 
(a) and YbCs (b) molecules. a is the Bohr length and is 
0.52918 A. The ground to excited state energy splitting in 
the separated limit was optimized to produce the Rb (Cs) 
5(6) s — > 5(6) p transition known experimentally. The arrows 
indicate the two-color photoassociation route to producing 
ultracold polar molecules in the ground X 2 E state. The inset 
shows the ground vibrational levels v — 0, 4, 8. 
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FIG. 2: (Color online) YbRb potential energy curves with 
spin-orbit effects included, ao is the Bohr length, 0.52918 A. 
The effective core potential ECP36SDF was modified to pro- 
duce the 5p 2 Pi/2 — + 5p 2 Ps/2 energy splitting of Rb. The ar- 
rows indicate the two-color photoassociation route. The inset 
shows the ground vibrational levels v — 0, 4, 8. 



TABLE II: Molecular parameters for YbRb including spin- 
orbital effects. Bond distance r e , dissociation energy (D), vi- 
bration constants (u> e and uj e Xe), rotation constants (B e and 
D e ), and vibration- rotation mixing constant a e are all deter- 
mined from the Morse potential parameter. r e is in atomic 
units while all others are in cm -1 . 



Molecule 


r e 


D 






B e 10~ 2 




a e 10 


YbRb X' z T, 1/2 


10.23 


193 


12.4 


0.20 


0.99 


25 


1.64 


YbRb a 2 n 1/2 


7.39 


3590 


70.9 


0.35 


1.89 


5.4 


1.00 


YbRb a 2 n 3/2 


7.52 


3120 


65.4 


0.34 


1.83 


5.7 


1.02 


YbRb fe 2 E 1/2 


8.58 


1026 


33.9 


0.28 


1.40 


9.6 


1.20 



orbital has a more pronounced presence in the second a- 
molecular orbital. There is some, albeit little, charge 
transfer from the alkali to Yb. At the potential minima, 
Rb transfers more charge to Yb than does Cs, counter to 
what is expected based on electro-negativity arguments. 
However, this can be understood in terms of the bond 
length in YbCs being larger than in YbRb. The excited 
a 2 II state of both YbRb and YbCs contains a far shorter 
bond length as well as a deeper well than the ground 
state X 2 S state. This is due to the larger polarizability 
in the 5(6) p state of Rb (Cs). 

In order to obtain a higher level of accuracy in de- 
scribing the excited a 2 IT state, we included the effects of 
spin-orbit coupling in the ab initio calculations. We then 
modified our long range potential to match the experi- 
mentally determined C§ coefficient determined by f]~3j ) . 
Fig. [II] gives the potentials with spin-orbit terms applied 
throughout. The spin-orbit terms of the ECP36SDF ef- 
fective core potential of the Stuttgart group were mod- 
ified so as to reproduce the 5p 2 Pi/ 2 — * 5p 2 P 3 / 2 energy 
splitting of the Rb atom to within a few cm -1 . The 
X 2 T,i/ 2 state is fairly unaffected by the presence of the 



spin-orbit terms; this is due to the large electronic energy 
separation of the ground and excited electronic states 
compared to the atomic spin-orbit constant in the alkali 
atoms. There are now 3 excited states since the 2 n state 
splits into two surfaces, one each for where the electronic 
spin lies parallel or anti-parallel to the electronic orbital 
angular momentum. In addition, at large separations of 
about 11 a there is an avoided crossing between the the 
° 2 ni/2 and the 6 2 E 1 / 2 states due to the spin-orbit in- 
teraction. This changes the outer wall behavior of both 
states. 

We again fit the ab initio points to a Morse potential in 
order to obtain the spectroscopic parameters describing 
the states of interest — see Table |TTJ There is very little 
change in the parameters for the X 2 T,i/ 2 state. The ex- 
cited a 2 n 1 / 2 state has a slightly smaller bond length but 
larger well depth while the the a 2 n 3 / 2 state has a longer 
bond length and smaller well depth compared to the a 2 IT 
state in the previous calculation. 



III. ALKALINE EARTH- YTTERBIUM 
MOLECULAR IONS 

In addition to the neutral molecules considered in the 
previous section, we can also consider their iso-electronic 
cation partners YbSr + and YbBa + . The molecules 
might be expected to form deeper wells and shorter 
bond lengths because the atomic cations Sr and Ba will 
strongly attract an electron. Also, these ions offer in- 
triguing possibilities for an EDM experiment in a sta- 
ble ion trap, as has been recently proposed!!?}. Another 
possibility is to arrange the molecular ions into an opti- 
cal lattice. In addition, replacing the alkali atoms with 
alkaline earths afford the possibility of isotopomers with 
zero nuclear spin on both atoms, thus greatly simplifying 
their spectroscopy. 

Weperformed a SCF-HF followed by an MCSCF 
[li . Il9l ] calculation on the neutral and ionic species of 
the alkaline earth- Yb separated atoms. Once again the 
Stuttgart basis set and ECPs ECP68MDF were used to 
describe the Yb 0] atom while the ECP36SDF and 
ECP54SDF[13] basis set and ECP were used to describe 
the Sr and Ba atoms respectively. CPPs were included as 
well. In the neutral molecule we accounted for the X 1 E, 
A 3 II, B 3 E, a 1 !! and b 1 ^ electronic states. In the ionic 
state we kept the three lowest 2 E states and the 2 U state. 
The final step in the calculations involved performing an 
MRCI0, El calculation on each spin multiplicity and 
group symmetry. The potential energy surfaces are pre- 
sented in Figs. |TTT] a and b. The upper (lower) set of 
curves in each panel are the ionic (neutral) molecules. 

A possible route to make the molecular ions is pre- 
sented in Figs. IHIl fa) and (b). By driving near the alka- 
line earth *S — > 1 P transition, one can photo-associate 
molecules into a high- lying vibrational level of the 1 II 
(or 1 E) state (solid arrow). From this state a photo- 
ionization pulse to the X 2 E state of the molecular ion 
should produce molecules in the electronic state of inter- 
est for an electron EDM experiment (long dashed poten- 
tial in the upper portion of each panel in Figs. IIIII (a) and 
(b)). 

A typical Ramsey experiment would measure the en- 
ergy difference between two levels with differing angular 
momentum projections M s referred to a laboratory mag- 
netic field. The measurement consists of the time- varying 
population of one of these states. To measure this, one 
could drive the molecule to the purely repulsive a 2 E state 
(short dashed arrow), which dissociates into the Sr (Ba) 
5(6) s 2 and Yb + 6s atomic states. If the dissociation laser 
were polarized to drive only one M s state, then efficient 
detection of the Yb + ions serves as a high signal-to-noise 
detection. 

These potential energy surfaces of the alkaline earth- 
Yb ionic molecules are markedly different from the neu- 
tral species alkali- Yb systems. First, the bond distance 
is reduced to around 8 a and the well depth increased 
by an order of magnitude. Physically, this is easy to un- 
derstand. The Sr and Ba ions (in the separated limit) 
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FIG. 3: (Color online) Potential energy surfaces for the YbSr 
(a) and YbBa (b) molecules in their neutral (lower set of 
curves) and ionic (upper set of curves) states, ao is the Bohr 
length, 0.52918 A. The ground neutral state to the ground 
ionic state energy splitting in the separated limit was op- 
timized to produce the ionization levels known experimen- 
tally. The arrows indicate the photo-association (solid ar- 
row), photo-ionization (medium dashed arrow) and proposed 
dissociation pulse (short dashed arrow) in the electron EDM 
search. 



are very willing to accept an electron, hence they have a 
tendency to take an electron from the neutral Yb atom. 
Thus, charge is efficiently transferred from the Yb atom 
to the alkaline earth ion and an ionic-type bond is formed, 
as compared to the more shallow van der Waals-type 
bond in YbRb and YbCs. This circumstance is also 
responsible for the larger body-fixed molecular electric 
dipole moment that in the neutrals. 

The second surface in the ion is the a 2 E state, which 
is repulsive above the equilibrium bond length of the 
A 2 E. The excited state surface in YbBa + is deeper than 
YbSr + , however the dissociation laser (DL) should still 
be able to drive the ground state to an unbound state 
of the system due to the rather large bond length of the 
a 2 E state in the molecular ions. The excited <z 2 E state 
is more closely akin to the ungerade states of homonu- 
clear molecules. In this instance, at the bond length of 
the A 2 E molecule, the configuration is a combination 
of two molecular orbitals, one of Ba(6s 1 )-Yb(6s 2 ) and 
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TABLE III: Molecular parameters for YbSr and YbBa. Bond 
distance r e , dissociation energy (D), vibration constants (w e 
and WeXe), rotation constants (B e and D e ), and vibration- 
rotation mixing constant a e are all determined from the 
Morse potential parameter. r e is in atomic units while all 
others are in cm -1 . 



Molecule 


r e 


D 


We 




B e IQ- 2 


D e IQ- 9 


Q e 10" 


YbSr X L H 


11.42 


94 


8.6 


0.19 


0.78 


26 


1.72 


YbSr a 1 II 


7.30 


8460 


94.0 


0.26 


1.98 


3.2 


0.63 


YbSr+ X 2 S 


7.94 


4230 


61.5 


0.22 


1.62 


4.5 


0.70 


YbBa X L H 


11.64 


112 


7.86 


0.14 


0.57 


12 


0.98 


YbBa a 1 II 


7.67 


6765 


73.4 


0.20 


1.32 


1.7 


0.41 


YbBa+ Y 2 S 


8.50 


2810 


42.8 


0.16 


1.08 


2.7 


0.47 



the other of Ba(6s 2 )-Yb(6s 1 ) hybridized atomic configu- 
rations. The excited state is the orthogonal combination. 
Whereas the ground state has an electronic distribution 
favorable to a deep bond the excited state does not. 

Table Mil gives the molecular parameters for a few 
states of interest in the YbSr and YbBa systems. Once 
again, we fit the tabulated points from the ab initio cal- 
culations to a Morse potential. The neutral molecules 
contain very shallow wells. This is attributed to each 
atom influencing the other and causing a slight mixing 
of atomic s- and p-orbitals. This is a van der Waals 
interaction which creates the ground neutral molecular 
states. However, the molecular ion states are created by 
having an electron from Yb have significant amplitude at 
the site of the alkaline earth atom. There is expected to 
be significant mixing of atomic s- and p-orbitals on both 
atoms. 



IV. EFFECTIVE ELECTRIC FIELD 
CALCULATION 

The effective electric field is defined by the energy shift 
incurred by the EDM Hamiltonian i/ c dm = — • £eS act- 
ing on a putative electric dipole moment d c . We have 
previously detailed a perturbative method for estimat- 
ing £ e g to within ~ 25% based on our non-relativistic 
molecular orbitals fl5l [28j, following the work in The 
basic idea is that we can separate the calculation into an 
analytic and numeric piece via (in atomic units) 

p (*mol|ffcdm|ttmol) 4(7 

where r ro i is a relativistic factor that depends on the 
atomic properties and is given by an analytic expression. 
(r = 1/2 is the electron projection on the molecular axis. 
Z is the number of protons in the heavy atom. The values 
of e s and e p represent the amount of s- and p-atomic or- 
bitals contained in the relevant molecular orbital, and are 
obtained from non-relativistic electronic structure calcu- 
lations. 



The estimate in the cases studied in [ll| [28| was for 
systems with one heavy atom and one light atom. There- 
fore, the dominant contribution to £ c g arises from the 
heavy atom. This is because the relativistic factor T rc i 
scales as Z 2 thus the effective electric field £ e s scales as 
Z a . In the examples we are considering here, both atoms 
contribute since they have similar Z. Therefore, it should 
be expected that the electron experiences an average ef- 
fective electric field from both atoms. What remains to 
be determined is whether the atoms are polarized paral- 
lel or anti-parallel to each other; this will be addressed 
below. 

Since r re i depends on the atomic properties of the atom 
considered, we break the two contributions up by expand- 
ing our molecular wave function as follows: 

|*mol> = £ s ,a|As) +e p , A |Ap) + 

£s,Yb|Ybs) +e PiYb |Ybp), (2) 

and verified the condition (^moil^moi) = 1 at all times in 
the calculation. The subscript A refers to either an alkali 
or alkaline earth atom. With this wave function we can 
calculate the quantity £ e g via 

<*mol|i?cdm|*mol) 



= £efi,A + £c«,Yb- (3) 

The only surviving terms are (A, s|Hedm|A, p) and 
(Yb, s|7J c dm|Yb, p) and so the contributions from each 
atom simply add. The terms are individually calculated 
via the methods described in [15j,|28|. Now, if the product 
£s,A£p,A is of the same (opposite) sign as e s ,Yb£p,Yb then 
the individual effective electric fields £ e ff,A and £ c g,Yb 
would add (subtract). Since each term has an uncer- 
tainty of 25% associated with it, and because these errors 
are not completely random, we report that these current 
estimates are 50% certain. 

One thing we discover in the ab initio calculations 
is that there are two large contributions to the overall 
molecular orbital, as discussed in Sec. [Til In the alka- 
lis, this mixing of molecular configurations is more pro- 
nounced than in the alkaline earths. These configura- 
tions are different zero angular momentum projections 
of atomic orbitals from each atom. We find that these 
two configurations are |^i) = \o-\aj) and |\& 2 ) = l cr 2 cr 3) 
for two different molecular orbitals o~i and 03. 

In the case of YbRb and YbCs 02 is primarily a Rb 
(Cs) 5s (6s) atomic orbital with appreciable mixing of 
Yb 6s and 6p z as well as Rb (Cs) 5p z (6p z ) atomic or- 
bitals. 0-3 has more of Rb (Cs) 5p z (6p z ) atomic charac- 
ter with larger admixtures of the other atomic orbitals. 
The YbSr + and YbBa + molecules also have similar dis- 
cernment in the molecular orbital 02- The 03 molecular 
orbital is composed of the alkaline and Yb p z atomic or- 
bitals with small admixtures of the respective s atomic 
orbitals. 

If we make the approximation that these two configu- 
rations are all that matter in the calculation of £ c g , we 
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TABLE IV: Effective electric field estimates for YbRb, YbCs, 
YbSr + , and YbBa + in GV/cm. Also presented are the molec- 
ular electric dipole moments (Debye) and the critical field for 
polarizing the molecules (kV/cm). 



Molecule 


£ e s 


See (A) 


fcff(Yb) 


dm 




YbRb X 2 E 


-0.70 


0.45 


-1.15 


0.21 


5.5 


YbCs A 2 E 


0.54 


1.42 


-0.88 


0.24 


3.5 


YbSr+ X 2 E 


-11.3 


10.6 


-21.9 


5.1 


0.38 


YbBa+ X 2 E 


1.2 


12.6 


-11.4 


5.1 


0.25 



can write the total molecular wave function as 



I ft 



= Ci|*i) +C 2 |*2 



(4) 



r 2 
< 2 



where c\ and c-i are the MRCI coefficients. Ideally, cf 
would be unity. However, since other configurations also 
contribute the sum c\ + C?, is approximately 0.95 in all 
the cases considered. 

The total s and p contribution from one of the atoms 
is given by the weighted sum (weighted by the MRCI 
coefficients) of the individual s- and p-atomic orbitals as 
follows: 



k 



3=1 



(5) 



In the above gj is a Gaussian centered on atom A hav- 
ing s-wave characteristics. The coefficients otj and fij 
describe the relative contributions of each Gaussian gj 
to the molecular orbital. c\ and c 2 are the same MRCI 
coefficients from the molecular orbital calculation. We 
can now define the e s a from atom A asfisl l28l 



(^A.sl^A,! 



(6) 



Similar definitions describe e p .A, £s,Yb, and e Pj Yb- The 
results of the calculation are presented in Table [TV] 

In Table IIVI we have separated out the contributions 
from both atoms individually. It is apparent that both 
atoms contribute to the size of the effective electric field, 
but with the opposite sign. We interpret this as both 
atoms polarizing in different directions with respect to 
the molecular axis. Indeed, this is what is observed when 
one looks at the electronic distribution of the unpaired 
electron. We also note that £ e ff of the YbRb and YbCs 
systems makes them an attractive alternative to experi- 
ments on the solitary alkali atoms. Even though £ e g is 
small compared to other X 2 Tj molecules, the expected 
long coherence times in an ultracold sample still makes 
them appealing. For the lone electron in Rb (Cs) atomic 
electron EDM experiments one would have to apply an 
electric field £ app of 27 MV/cm (4.5 MV/cm) — using 
the enhancement factors in Q — to achieve the same 
£ c ff- This is to be contrasted with the proposed alkali- 
Yb molecular electron EDM experiment where one needs 
£app ~10-25 kV/cm. The values of £ app in the molecular 



cases are based on multiplying £ po i by 3-5 in order to be 
in the linear Stark regime. £ po \ is defined as the electric 
field required to equalize the Stark energy and the energy 
splitting between rotational levels in the molecule. £ po i 
is tabulated in Table [TVl 

We have used the following convention for the overall 
sign of £ Q g: we define the molecular axis n as pointing 
from Yb to A, where A is either an alkali or alkaline earth 
atom. This is because we choose the axis to point from 
more negative charge to more positive. In the alkali- Yb 
molecules, Yb has negative charge while in the alkaline 
earth- Yb molecules it has less positive charge. A positive 
value of £ c ff means that it lies against the direction of n, 
i.e. points from positive charge toward negative charge. 
Positive feff, in turn, means that the atomic electron is 
polarized so that it is displaced in the direction — h. 

These molecules offer alternatives to performing the 
experiments on Rb and Cs or Sr + and Ba + alone. In the 
atoms, the size of £ s is proportional to the applied field 
£app by an enhancement factorQ. In molecules, this is 
not the case provided they are fully polarized. However, 
polarizing these molecules is not as easy to do as in the 
case of 3 Ai molecules discussed in fl5i |28| . Yet, the 
£ e g that can be attained is larger than the £ e ff that can 
be attained in the atomic systems for a similar £ app 0]. 
Therefore, with the recent advances in producing cold 
molecules, these systems offer very intriguing and viable 
alternative routes to probe for the electron EDM. 

The idea of creating ultracold molecules from ultracold 
atoms, along with the anticipated benefit of molecules 
to electron EDM searches, to explore a variety of laser 
cooled atoms paired with Yb for an electron EDM exper- 
iment is the focus of this paper. We chose heavier pairing 
partners in the hope of reducing the rotational constant 
of the molecules, thereby making them easier to polar- 
ize in modest applied electric fields. However, as one in- 
creases the mass of the pairing partner, one also increases 
the contribution it gives to £ g and this contribution is 
always opposite to Yb's contribution, thereby reducing 
the effect. Therefore, there should be some trade off be- 
tween reduced rotational constant and overall £ e ff- Fu- 
ture work would include exploring systems such as YbK, 
YbNa, YbCa + , and YbMg + . In addition, the same ideas 
would work for Hg-paired partners. Work on HgH[29[ has 
already been performed in regards to an electron EDM 
experiment. Systems such as HgLi, HgNa, etc. should 
also give similar results as the Yb-paired systems with 
the added benefit that Hg is heavier than Yb, thereby 
hopefully producing a larger £ e g. 



V. CONCLUSIONS 

We have considered a particular class of 2 S molecules 
that offer some advantages and drawbacks to current 
polar molecule electron EDM searches. This type of 
molecule is formed from ultracold atoms so that the 
ground state molecules are also ultracold, leading to the 
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possibility of longer coherence times. We have consid- 
ered YbRb (-0.7 GV/cm), YbCs (+0.54 GV/cm), YbSr+ 
(-11.3 GV/cm) and YbBa+ (+1.2 GV/cm). 

YbRb and YbCs can be polarized in fairly modest 
fields and thereby are intriguing inasmuch as they of- 
fer some incentive over searching for the electron EDM 
in the alkali atom alone. The molecular ions YbSr+ and 
YbBa + have many useful attributes, However, the field 
required to polarize these molecules is rather large for 
an ion trap. Thus, the use of an optical lattice may be 
more practical for holding these ions for a long period 



of time without having to rotate an £ app on the order of 
1 kV/cm. 
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